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Исследование поведения нестационарного 
динамического объекта в бортовых ИС 
новых поколений” 


Рассмотрен подход, основанный на реализации принципа конкуренции в рамках концепции мягких 
вычислений. Методологической основой для построения алгоритмических процедур принятия 
решений в нечетких условиях является концепция оптимальности решения и имеющихся ресурсов 
(нечеткая цель и ограничения). 


Введение 


Нестационарные системы — достаточно распространенный класс динамических 
систем. Изучение поведения таких систем представляет большой научный и практи- 
ческий интерес при контроле динамики сложных объектов. В настоящей работе 
развивается подход к синтезу нейросетевых моделей систем, характеристики которых 
изменяются с течением времени. Оценка параметров и определение структуры модели 
осуществляются на основе разработанного алгоритма структурного и параметрического 
синтеза нейронной сети. 

Эффективность процедур обработки информации и принятия решений в значи- 
тельной степени определяется качеством моделей, используемых как для описания 
объектов и закономерностей, так и для прогнозирования будущего состояния систем. 
Во многих прикладных задачах возникает необходимость в обработке нестационарных 
последовательностей и формировании модели систем, характеристики которых из- 
меняются с течением времени. Одной из таких задач, связанных с использованием 
нейросетевых технологий, является контроль характеристик морского динамического 
объекта (ДО) в условиях непрерывного изменения параметров объекта и внешней 
среды [1-26]. 

Типичные картины возникновения и развитие нестационарных процессов мож- 
но наблюдать при изучении процессов взаимодействия аварийного ДО с внешней 
средой при контроле непотопляемости судна в бортовой интеллектуальной системе 
(ИС) «Непотопляемость». Поступление воды в отсеки ДО и ее фильтрация в сосед- 
ние отсеки при получении пробоины — типичные картины возникновения и развития 
нестационарных процессов в сложных динамических системах. Основная задача, ре- 
шаемая в процессе функционирования ИС, состоит в выработке практических реко- 
мендаций по прогнозу развития аварийных ситуаций в условиях неопределенности и 


* Исследования проведены при поддержке Федерального агентства по образованию по проекту 
«Интеллектуальная система навигации и управления морским динамическим объектом в экстремальных 
условиях эксплуатации», в рамках реализации ФЦП «Научные и научно-педагогические кадры иннова- 
ционной России» на 2009 — 2013 годы. 
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неполноты исходной информации. Для реализации этой задачи авторами предложен 
подход к синтезу нейросетевых моделей нестационарных систем, включающий проце- 
дуры структурного и параметрического синтеза ИНС, а также параметров распределений 
на основе специально разработанного алгоритма. При использовании синтезированных 
моделей для обработки последовательно поступающих данных измерений выполняется 
рекуррентная модификация параметров нейронной сети. Разработанные процедуры могут 
быть использованы при решении ряда задач идентификации системы (определения теку- 
щей ситуации) и краткосрочного прогнозирования. 


Математическая модель, описывающая физические 
картины поведения аварийного динамического объекта 


Математическая модель, описывающая взаимодействие аварийного ДО с внешней 
средой, представлена системой дифференциальных уравнений бортовой качки судна 
с затопленным отсеками [4], [14], [20]: 


пер» Св: =Р“(В)} (1) 


п=1 


(т + из.) пе + + Ат 


(И +и.)9"+ 4.909 |9 ]+ Б"Г (9) + ши + Ато + 

+рУ (а, +8С,в,)=М" +" (1 (Е -21)-ЛА" 0; и 
ево, ав =) а 
+5С, (9+4")+ С, (т. +п"")} =0, (3) 


(п — В) 
где с = (. -п’) — относительное поперечно-горизонтальное перемещение судна; 


9 = (@ — №) — относительный угол крена; 2’ — аппликата центра тяжести затопленного 


объема; знаком «*» отмечены величины, претерпевшие вариацию в результате затопле- 
НИЯ. 

Интерпретация нештатных ситуаций дается на основе представлений в виде 
временных кривых и траекторий изображающей точки на фазовой плоскости. 
Анализ поведения аварийного ДО на волнении позволяет выделить три характерные 
ситуации, определяющие развитие и стабилизацию нестационарных процессов бор- 
товой качки аварийного ДО: 

— развитие и стабилизация крена (после проведения мероприятий по борьбе за 
живучесть) в процессе взаимодействия ДО с внешней средой на нерегулярном 
волнении ограниченной интенсивности; 

— переход ДО из одного стабильного состояния в другое (смена аттрактора) 
после завершения процесса затопления отсеков и прекращения фильтрации воды в 
соседние отсеки; 

— непрерывное нарастание динамического крена при обвальном затоплении 
отсеков аварийного судна, находящегося под воздействием нерегулярного волнения. 

На рис. | - 4, 6, 7 представлены экстремальные ситуации, характеризующие 
указанные выше режимы колебаний. Первая из этих ситуаций (рис. 1) возникает при 
симметричном затоплении отсеков и фильтрации воды из смежных отсеков асим- 
метричного затопления, вследствие чего ДО приобретает крен (в рассмотренном сце- 
нарии развития ситуации 6 = 12°) и находится под воздействием внешних возму- 
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щений (ветер, нерегулярное волнение). Кренящий момент от фильтрации воды в 
асимметричные отсеки изменяется монотонно до угла 6 = 12°, затем стабилизируется, а 
после проведения мероприятий по борьбе за живучесть монотонно убывает и 
система возвращается в исходное состояние: 


М (1) =20(0<12°1<1), М (1) = сопяЕ (0 =12°,1 в [1,1,]), М (1) =Ю01 = [Ь,4]. 


В результате проведения мероприятий по борьбе за живучесть (откачка посту- 
пившей воды, прекращение фильтрации воды из смежных отсеков, устранение крена 
до 0 = 0°) ДО возвращается в устойчивый (близкий к исходному) режим колебаний 
на нерегулярном волнении. 

[№4 моделирование ЕЕ" х^— 


Выход Исходные данные Вычисления Дополнительно 


Исходные данные [График угла бортовой качки || Фазовая траектория | Плечо восстанавливающего момента | Г_4 | > 


Теа 
Че 1 


1 200 


Рисунок 1 — Временные кривые, характеризующие развитие и стабилизацию крена 
при асимметричном затоплении отсеков ДО на нерегулярном волнении (светлая 
кривая описывает эволюцию угловой скорости при бортовой качке) 


Фазовый портрет рассматриваемой динамической системы представлен на рис. 2, 
из которого следует, что формирование аттракторного множества в этом случае 
представляет собой типичный случай одностабильной динамической системы, под- 
верженной нестационарным случайным возмущениям, выводящим ДО из стабиль- 
ного состояния на определенном отрезке времени, и вновь возвращаемой в исходное 
состояние после снятия возмущения. 

ЧЕ Моделирование ее" х— 


Выход Исходные данные Вычисления Дополнительно 


Исходные данные | График угла бортовой качки | Фазовая траектория || Плечо восстанавливающего момента | Г_4 |» 


Оастаау обаавобёу 


Рисунок 2 — Фазовый портрет, характеризующий поведение динамической системы в 
процессе развития и стабилизации крена 


Вторая ситуация (рис. 3) связана с переходом динамической системы от одного 
стабильного состояния равновесия (0 = 0) к другому 0 = 61, более неблагоприятному 
стабильному состоянию (смена аттрактора). Физические особенности такого перехода 
объясняются возникновением крена ДО вследствие асимметричного затопления отсеков. 


512 «Искусственный интеллект» 372010 


Исследование поведения нестационарного динамического объекта... БН 


Как видно из рис. 4, при крене ДО 090 =0, система представляет собой последова- 
тельность нерегулярных колебаний около угла 9! = 17° под влиянием различных (внут- 
ренних и внешних) факторов и, в первую очередь, нестационарного ветра и нерегу- 
лярного волнения. 


8 Моделирование [ею] 
[ Выход _ Исходные данные Вычисления Дополнительно 
| 

`Исходные данные [График угла бортовой качки | Фазовая траектория | Плечо восстанавливающего момента | Г_4 | > 


овой качки и Производная изменения угла бортовой качки 


Теа 
ава 


Рисунок 3 — Смена стабильных состояний (смена аттрактора) динамической системы 


Фазовый портрет динамической системы, представленной на рис. 3, характерен для 
многих нелинейных систем, описываемых с помощью математической модели (1). 
Эволюция фазовых траекторий аварийного ДО при асимметричном затоплении отсеков 
приведена на рис. 4. 


[= [151888389] 


9 Моделирование 
| Выход Исходныеденные Вычисления Дополнительно 


Исходные данные | График угла бортовой качки |Фазовая траектория | Плечо восстанавливающего момента | Г_4 |? 


‘бастаау обаавотову, 


0,069  -0,009 — 0,051 0,112 0.172 — 0,232 0,292 0,352 0,412 0,472 0,532 


Рисунок 4 — Фазовый портрет динамической системы, характеризующей смену 
стабильных состояний ДО на нерегулярном волнении 


Из этого рисунка следует, что формирование аттракторного множества в этом 
случае представляет собой типичный случай бистабильной динамической системы, 
характеризуемой нестационарным поведением. Случайные возмущения (нерегулярные 
колебания жидкости в затопленных отсеках, нестационарный ветер и нерегулярное вол- 
нение), выводят ДО из стабильного состояния при 6 = 0 в новое состояние при 0 = 01. 

Следует заметить, что представленный на рис. 4 аттрактор отличается от 
классического аттрактора бистабильной системы (рис. 5 В), формирование фазовых 
траекторий в которой вызвано особенностями динамической системы (третий клас- 
сический случай затопления). Эти особенности можно объяснить тем, что исследуемый 
ДО является неустойчивым в прямом положении. Он может наклоняться на правый (ВВ) 
или левый (ГВ) борт, (причем @ьв = 0:в) даже при отсутствии внешнего возмущения 
(случай отрицательной метацентрической высоты, вызванной значительной свободной 
поверхностью жидкости в затопленных отсеках). Одним из факторов, воздействие 
которых может привести к возникновению ситуации, представленной на рис. 5 В), 
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является смещение груза в трюмах, другим — потеря начальной остойчивости при 
затоплении отсеков. Именно последняя ситуация является типичным примером ди- 
намики бистабильной системы в классическом ее понимании (рис. 5 В). 

0,8 


Рисунок 5 — Фазовый портрет трехстабильной (А) и бистабильной (В) динамической 
системы при анализе экстремальных ситуаций 


В этом случае в динамической системе имеются два аттрактора, и нерегулярные 
колебания происходят около равновесных, симметрично расположенных углов крена 
на правый или левый борт. Еще более сложный случай — поведение трехстабильной 
системы (рис. 5 А), которая характеризуется наличием трех аттракторных множеств. 
Возникновение такой сложной колебательной системы возможно в случае, когда ДО 
имеет развитые надстройки, водонепроницаемость которых обеспечена в процессе 
эксплуатации [11]. Переход от ситуации (рис. 3) к ситуации (рис. 5 В) возможен при 
условии потери начальной остойчивости в условиях симметричного затопления 
больших отсеков, расположенных от борта до борта (достаточно типичный случай 
для судов небольшого водоизмещения). 

Третья ситуация (рис. 6) характеризуется непрерывным нарастанием крена в 
результате поступления больших масс воды внутрь корпуса ДО. Поведение такой 
динамической системы представляет типичный случай нестационарных колебаний и 
связано с непрерывным изменением динамики объекта и внешней среды. 

Это — наиболее тяжелый случай реализации динамики аварийного объекта на 
волнении. Из-за поступления больших масс воды в условиях асимметричного за- 
топления ДО теряет остойчивость и опрокидывается под воздействием внешних 
возмущений. 


[148 Моделирование ЕЕ х—) 


Выход Исходные данные Вычисления Дополнительно 


Исходные данные График угла бортовой качки | Фазовая траектория | Плечо восстанавливающего момента | Г_4 | > 


Рисунок 6 — Непрерывное нарастание крена при обвальном затоплении отсеков 
аварийного ДО на волнении 
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Фазовый портрет нестационарной динамической системы, представленной на 
рис. 6, изображен на рис. 7 и характеризует непрерывное нарастание крена, сопро- 
вождающееся нерегулярными колебаниями ДО. В процессе эволюции рассматриваемая 
система либо окажется устойчивой (перейдет в новый аттрактор) при постепенном 
уменьшении интенсивности затопления и внешних возмущений (как это видно из 
рис. 3 и 4), либо произойдет потеря устойчивости колебательного движения (опроки- 
дывание, рис. 6 и 7, когда колебательный режим реализуется около угла 0 = 180°). 
=. 88:38) 


18 Моделирование 


Выход Исходные данные Вычисления Дополнительно 


Исходные данные ] График угла бортовой качки [Фазовая траектория |] плечо восстанавливающего момента ] г4|› 


А 


-0,062 0,050 0,161 0,273. 0,385 0,497 0,609 0,721 0,833 0,945 1,057 


Рисунок 7 — Фазовый портрет нелинейной нестационарной динамической системы, 
характеризующий потерю устойчивости нерегулярных колебаний в условиях 
непрерывного изменения динамики объекта и внешней среды 


Геометрическая интерпретация 


Рассмотрим принципиально иной подход к моделированию нелинейной нестацио- 
нарной динамической системы для рассмотренных экстремальных ситуаций. Используем 
методы теории катастроф [23] и представим поведение системы с помощью ката- 
строфы сборки, позволяющей получить наиболее полное представление о топологии 
геометрического пространства в процессе развития эволюции системы. 

Обсудим сначала первые две стратегии поведения нелинейной нестационарной 
системы с позиций теории катастроф. На рис. 8 — 9 представлены топологические 
картины взаимодействия аварийного судна с внешней средой. 


Рисунок 8 — Топологические картины поведения нелинейной нестационарной 
системы в экстремальной ситуации: А) — первый тип затопления; 
В) - второй тип затопления 
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Здесь изображены картины модифицированной катастрофы сборки применительно 
к исследуемой физической системе, показаны характерные точки, кривые и области 
изменения исследуемых характеристик в плоскости У, для экстремальной ситуации, 
характеризующей динамику аварийного судна на волнении: пло — исходное положение 
метацентра; 11(0) — метацентрическая эволюта (многообразие катастроф); Си — исход- 
ное положение ЦТ. Заштрихованные области С7(0) характеризуют конечные множе- 
ства, отображающие динамическую среду в конечных стадиях развития аварийной 
ситуации. 

Для аналитической интерпретации ситуации в правом верхнем углу представлена 
диаграмма изменения восстанавливающего момента М(6). Цифрами 1, 2, 3 на рис. 8 А) 
отмечены кривые М(0) для положений ЦТ в точках Со, С: и С› соответственно. 
Точка Си является центром величины (ЦВ) судна при угле крена 0 =0. Кривая С(6) 
представляет собой множество, отображающее структурные изменения в форме под- 
водной части корпуса. С помощью этой кривой описывается пространственное пере- 
мещение ЦВ в зависимости от 0. Радиус кривизны кривой С(0) определяет величину 
метацентрического радиуса г(6), величина которого при 6 = 0 характеризует начальное 
положение метацентра то. Кривая метацентров образует метацентрическую эволюту — 
бифуркационное множество В(0) модифицированной модели сборки. В силу сим- 
метрии кривая С(0) всегда будет иметь при 9 =0 либо локальный минимум, либо 
локальный максимум, а кривая метацентров — особенность, в данном случае особен- 
ность стандартной или двойственной сборки. 

Интерпретация ситуации на рис. 8 А) дает следующие результаты. В момент 
появления пробоины система скачком перемещается в новое состояние. Возникновение 
скачка (катастрофы) обусловлено резким ухудшением начальной метацентрической 
высоты вследствие принятия жидкого груза и образования большой площади 
свободной поверхности жидкости симметрично затопленного отсека. ЦТ системы 
практически мгновенно перемещается из точки Су в точку С1. При этом значительно 
уменьшается область С7/(0) и восстанавливающий момент М(6) и судно оказывается 
в критическом состоянии, поскольку интенсивные внешние возмущения могут при- 
вести к опрокидыванию судна. После проведения мероприятий по стабилизации 
ситуации система возвращается в область с ЦТ в точке С2. При этом остается незна- 
чительный неуравновешенный крен 091. Нестационарность исследуемых динамических 
процессов в рассматриваемой ситуации проявляется в непрерывном изменении 
математического ожидания и дисперсии в процессе эволюции системы. 

Топологическая картина на рис. 8 может быть дополнена другими точками, 
позволяющими рассматривать различные состояния судна и определяющими экви- 
валентные положения ЦТ для опрокидывающего момента, углов крена, соответству- 
ющих начальному (во время приложения кренящего момента) и конечному (дина- 
мическому) наклонению судна, а также эквивалентных положений ЦТ твердого и 
жидкого груза во время погрузки (разгрузки) судна. 

Случай на рис. 8 В) характеризует асимметричное затопление. Здесь происходит 
пространственное перемещение ЦТ из точки Оо в точку С!. Скачок возникает при 
достаточно интенсивном затоплении. Если удается стабилизировать ситуацию, то ЦТ 
сохраняется в точке С1, а бифуркационное множество В определяется положением 
эволюты с точкой сборки Р. Аналитическая интерпретация ситуации может быть 
проведена с использованием диаграммы остойчивости в правом верхнем углу экрана 
для положения ЦТ в точке (2. Из этих данных следует, что в рассматриваемом 
случае затопления судно плавает с креном 0! в условиях пониженной остойчивости. 
Это требует проведения неотложных мероприятий по поддержанию живучести судна. 
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Рисунок 9 — Топологическая картина развития экстремальной ситуации, 
описывающей поведение нелинейной нестационарной системы в условиях 
непрерывного изменения динамики объекта и внешней среды 


Ситуация на рис. 9 еще более тяжелая. Из-за непрерывного поступления внутрь 
корпуса больших масс воды происходит непрерывное изменение динамики объекта. 
Причем, точка С! здесь соответствует более высокому положению ЦТ по сравнению 
с ситуацией на рис. 8, а диаграмма остойчивости характеризуется кривой М(0), свиде- 
тельствующей о чрезвычайно низкой сопротивляемости к восприятию кренящих на- 
грузок. В результате судно оказывается в критическом состоянии, поскольку доста- 
точно незначительного кренящего момента, чтобы вызвать опрокидывание. В случае 
фильтрации воды в соседние отсеки положение аварийного судна еще более ухуд- 
шится и потребует срочного принятия решений по стабилизации ситуации. 

Если продолжается обвальное затопление (рис. 9), то система непрерывно 
изменяет свое состояние в процессе эволюции и совершает движение в направлении 
возникновения новой катастрофы — опрокидыванию под воздействием внешних 
возмущений в условиях резкого снижения практически всех элементов, определяющих 
катастрофу сборки. 


Заключение 


Анализ использования нейросетевых технологий при исследовании поведения 
аварийного ДО как нестационарной динамической системы в условиях нерегулярного 
волнения показал, что процесс взаимодействия объекта с динамической внешней 
средой характеризуется рядом особенностей, присущих сложным нелинейным систе- 
мам, функционирующим в условиях неопределенности и неполноты исходной ин- 
формации. Трудности, связанные с анализом таких систем, в значительной степени 
могут быть преодолены за счет применения нейросетевых технологий. Рациональное 
использование ИНС при исследовании нестационарных динамических систем позво- 
ляет обеспечить гибкость и способность адаптироваться к изменяющимся внешним 
условиям. При этом сохраняется устойчиво высокое качество работы ИНС и реализу- 
ется недосягаемый ранее уровень вычислительной мощности. 

Проведенные вычислительные эксперименты свидетельствуют о работоспособ- 
ности и эффективности реализации нейросетевых моделей в сложных задачах анализа 
и интерпретации экстремальных ситуаций в условиях непрерывного изменения 
динамики объекта и внешней среды. Параллелизм и возможность нелинейного преобра- 
зования информации в реализуемых конфигурациях нейросетевых структур открывают 
перспективы практической реализации ИНС в задачах анализа и интерпретации 
данных динамических измерений в бортовых ИС новых поколений. 
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о 


Ю.Л. Нечаев, О.М. Петров 

Дослдження поведшки нестацюнарного динам1чного об’екта в бортових ГС нових поколнь 
Розглянуто шджд, побудований на реалзацй принципу конкуренцй у рамках концепций м’яких обчислень. 
Методолоичною основою для побудови алгоритм!чних процедур прийняття рипень у нечтких умовах е 
концепщя оптимальности! рипення та наявних ресурслв (нечтка мета та обмеження). 
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